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We study atom-atom correlations and relative number squeezing in the dissociation of a Bose- 
Einstein condensate (BEC) of molecular dimers made of either bosonic or fermionic atom pairs. Our 
treatment addresses the role of the spatial inhomogeneity of the molecular BEC on the strength of 
correlations in the short time limit. We obtain explicit analytic results for the density-density 
correlation functions in momentum space, and show that the correlation widths and the degree of 
relative number squeezing are determined merely by the shape of the molecular condensate. 
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Dissociation of a Bose-Einstein condensate (BEC) of 
molecular dimers [l| into pair-correlated atoms repre- 
sents the matter-wave analog of two-photon parametric 
down-conversion. The latter process has been of cru- 
cial importance to the development of quantum optics. 
Owing to this analogy, molecular dissociation currently 
represents one of the "workhorses" of the new field of 
quantum-atom optics (3, [3, 0, H, [3, 0, @, [3, E3] and offers 
promising opportunities for the generation of strongly 
correlated atomic ensembles and fundamental tests of 
quantum mechanics with macroscopic numbers of mas- 
sive particles. Examples include the demonstration of 
the Einstein-Podolski-Rosen paradox and violation of a 
classical Cauchy-Schwartz inequality [3, [3, EH ■ A closely 
related process is atomic four-wave mixing in a collision 
of two BECs [13, E3], which produces a spherical halo of 
spontaneously scattered atoms 3, 15, la, 12, 13, 13, 23] 
with correlations very similar to those in dissociation. 

A recently emer ged discussion topic - following the ex- 
periments of Ref. [13j on BEC collisions - is the under- 
standing of the width and the strength of the observed 
correlations, as well as the prospects of detecting rela- 
tive number squeezing in the halo of the s-wave scat- 
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20] . The same questions are relevant 
to atom-atom correlations in molecular dissociation and 
have not been fully addressed so far. 

In this paper, we study atom-atom correlations and 
relative number squeezing in the dissociation of a molec- 
ular BEC in the short time limit. Our analysis applies 
to molecules that may consist of pairs of either bosonic 
or fermionic atoms, and takes into account the spatial 
inhomogeneity of the molecular BEC. It has been ar- 
gued before and shown in the bosonic case @, [3] that 
the treatment of spatial inhomogeneity is crucial for ob- 
taining quantitatively correct results for atom-atom cor- 
relations. In the fermionic case, the treatments of dissoci- 
ation have been so far restricted only to uniform systems 
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[22|. In all cases, however, the techniques are 
numerical and do not give the transparency of analytic 
understanding, in contrast to the results obtained here. 



The effective quantum field theory Hamiltonian de- 
scribing our system in the undepleted molecular field ap- 
proximation is given, in a rotating frame, by [23] 
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The operators ^>i^{x, t) describe the atoms in two differ- 
ent spin states, which can be either bosonic or fermionic, 
and we assume that they have the same mass. The effec- 
tive coupling g(x) is defined as g(x) = \\/ Po(x), where 
X is the atom-molecule coupling (see [22| for details) and 
Po(x) is the initial density of the molecular BEC in a 
harmonic trap. For computational simplicity we start 
by treating a one-dimensional (ID) system; the analytic 
results will later be generalized to three dimensions (3D). 

The key difference between the present and previ- 
ous (uniform) treatments of dissociation in the unde- 
pleted molecular approximation 0, [3] is that we retain 
the spatial dependency of the molecular BEC: the effec- 
tive coupling g(x) absorbs the molecular field, which is 
treated classically via the coherent mean-field amplitude 
*&o(x) = \J po (x) . The undepleted molecular approxi- 
mation is valid only for short dissociation times, during 
which the converted fraction of molecules does not exceed 
~ 10% [zlH^]- Accordingly, the coupling g{x) can be kept 
constant in time, although the evolution of the atomic 
field is taking place in free space. In this regime, the dis- 
sociation typically creates low-density atomic clouds for 
which the s-wave scattering interactions are a negligible 
effect too 0]. 

The detuning A in Eq. ([1]) corresponds to the energy 
mismatch 2hA between the free two-atom state in the 
dissociation threshold and the bound molecular state. 
Molecules that are unstable against spontaneous disso- 
ciation correspond to A < 0, with 2?i|A| being the total 
dissociation energy that is converted into kinetic energy 
of atom pairs primarily populating the resonant momenta 
around ±/co, with ko = y/2m\A\/h. 
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Writing down the Heisenberg equations of motion 
for the field operators and converting to Fourier space 
tyj(x,t) = J dkaj(k,t)exp(ikx)/^/2ir, we arrive at the 
following coupled equations for the operators a,j(k,t): 

dcii(k, t) 
dt 

da\(k, t) 
dt ~ 



-7=g(q + k) a\ (q,t), 

V Z7T 



iA k al(k,t) 



dq 



ff(g-A)oi(-«,t). (2) 



Here and hereafter the + (— ) [in general, upper (lower)] 
sign stands for bosonic (fermionic) statistics of the atoms, 
g(k) = J dxe~ tkx g{x) / ^/2~n is the Fourier transform of 
the effective coupling g(x), and A k = hk 2 / (2mi) + A. 

Equations ^ can be solved numerically using stan- 
dard methods of linear operator algebra. One can show 
that, for vacuum initial conditions, the only nonzero 
second-order moments are the normal and anomalous 
atomic densities, rij(k,k',t) = (a](fc,t)a.,'(fc / ,t)) and 
mi2(k,k',t) = (ai(k,t)a 2 (k' ,t)) . From Eqs. <[2j) we can 
see that the finite width of g(k) - due to the inhomo- 
geneity of the source - implies that ai(k) couples not 
only to the partner spin component at exactly opposite 
momentum a\{—k) (as is the case in the homogeneous 
system), but also to a range of momenta around —k, 
within — k ± 8k. The spread in 5k determines the width 
of atom-atom correlations and is ultimately related to 
the width of the molecular BEC. 

We now turn to the quantitative analysis of atom-atom 
correlations expected to be present between the differ- 
ent spin-state atoms with equal but opposite momenta 
due to momentum conservation, and between the same 
spin-state atoms in the collinear direction due to quan- 
tum statistical effects. We quantify these correlations via 
Glauber's second-order correlation function, 

(2) _ (al(k,t)a](k',t)a j (k',t)a i (k,t)) 

9 » [ ' ' ' ~ ni (k,t) nj (k',t) ■ (Jj 

The normalization with respect to the product of den- 
sities rii(k,t) and rij(k',t) [with rij(k,t) = rij(k, k,t)] 

(2) 

ensures that g\- (k,k l ,£) = 1 for uncorrelated states. 

(2) 

Due to obvious symmetry considerations, g{ 2 (k, & 3 1) — 
g$> (k, k', t) and g$ (k, k',t) = g$ (k, k',t). 

Since the effective Hamiltonian corresponding to 
Eqs. ([2]) is quadratic in the field operators, we can ap- 
ply Wick's theorem to factorize the fourth-order mo- 
ment in Eq. ©. Noting that (a{(k, t)a 2 {k' , t)) = 
(a,j(k,t)a,j(k' ,<)) = in the present model, we obtain 

g^ (k, k', t) = l + \m 12 (k, k', t)\ 2 /[ ni (k, t)n 2 (k', t)}, (4) 

gf){k,k',t) = l±\n j (k,k',t)\ 2 /[n j (k,t)n j (k',t)}. (5) 

Before presenting the results based on numerical so- 
lutions of Eqs. ([2]), we now develop simple analytic ap- 
proaches that give approximate predictions for these ob- 
servables, valid for short times. More specifically, we 



treat the short time dynamics of dissociation via the Tay- 
lor expansion in time, up to terms of order t 2 [26], 



a,j(k, t) — a,j(k, 0)- 



ddj(k, t) 
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valid for t <C to, where to — l/X\/po(0) is the time scale. 
Using the rhs of Eqs. ([2]), this gives, up to the lowest- 
order terms, rij(k,k',t) ~ t 2 J dqg*(q + k)g(q + k')j2-K 
and \mi2(k, k' , t)\ ~ t\g(k + k')\/\f2n, or equivalently 

nj (k, k', t)~t 2 [ dxe-^ k - k "^ x [g{x)} 2 /2n, (6) 



dxe- l{k+k ' )x g{x)/2n . (7) 



\m 12 (k,k',t)\ ~ t 



These results show that the width of the collinear (CL) 
correlation between the same-spin atoms with nearly the 
same momenta, Eq. J5]) , will be determined by the square 
of the Fourier transform of the square of the effective cou- 
pling g{x). On the other hand, the width of the back- 
to-back (BB) correlation, Eq. Q, between the different 
spin-state atoms with nearly opposite momenta will be 
determined by the square of the Fourier transform of 
g(x). Therefore, the CL correlation is generally broader 
than the BB correlation. These conclusions are true for 
any shape of the source and apply to both bosonic and 
fermionic statistics in the short time limit. 

Thomas-Fermi (TF) parabolic density profile. - We 
now give explicit analytic results for the case of a TF in- 
verted parabola for the molecular BEC density, po{x) — 
po(l ~ x 2 / R^-p) (\x\ < Rtf), in which case g{x) = 
X v / po(l — x 2 /R 2 ?F ) 1 / 2 . U sing the integral representation 
of Bessel functions J v (z) (ij, Eqs. © and (0 yield 

^k^^^^^Zp^, ( 8 ) 
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Since J v {z) ~ (z /2) u /T{v + 1) for z«l, the atomic mo- 
mentum distribution rij(k, t) and the diagonal anomalous 
density mi 2 {k,—k,t) are rij(k,t) ~ 2t 2 x 2 PoRtf /3ir and 
\mi 2 (k,—k,t) \ ~ txyfpoRTF / 4. Despite the fact that the 
atomic momentum distribution in the lowest order in t is 
uniform, the momentum cutoff k max [23| - which must be 
assumed when using a (5-function interaction in Eq. ([1]) - 
prevents the total atom number from diverging. 

Substituting Eqs. ® and © into Eqs. (g]) and (0), we 
obtain the following explicit results for the atom-atom 
correlations, valid for t <C to: 
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Figure 1: Back-to-back (a) and collinear (b) correlation 

(2) 

glj (k,ko,t) as a function of k at t/to = 0.5. The dimension- 
less detuning is 8 — At = —9, where to = 1/Xv Po(0) ~ 5 ms 
is the time scale; the Thomas-Fermi radius of the molecular 
BEC is Rtf = 250 fim; for other parameters, see Ref. [24|. 



(2) 

The pair correlations <?■• '{k,ko,t = 0.5tn), where the 
momentum of one of the atomic components is fixed to 
ko, while the partner momentum k is varied, is plotted 
in Figs. UK a) andQJb). The dashed lines are the ana- 
lytic results of Eqs. fTTjj) and (fTTj) . whereas the solid lines 
are the numerical results from Eqs. (J2j) - In Fig.QJb) the 
dashed lines are almost indistinguishable from the respec- 
tive solid lines, even though the ratio t/to — 0.5 is not 
very small. For earlier times, the agreement between the 
analytic and numerical results is even better. 

( 2) 

In the case of different spin states, g\ 2 J (k,ko,t), with 
either bosonic or fermionic atoms, we see a strong BB 
correlation signal between atom pairs with equal but op- 
posite momenta, centered at k = — k . The same-spin 
CL correlation function ■ (k, ko,t), on the other hand, 
shows the Hanbury Brown and Twiss bunching peak for 
bosons and an antibunching dip for fermions due to Pauli 
blocking [iij], both centered at k = k - 

In Fig.^a), we plot the widths of the BB and CL cor- 
relations as a function of time. For simplicity, we define 
them as the half-width at half-maximum. The widths 
in the bosonic and fermionic cases (with the solid lines 
corresponding to the numerical results) have universal 
asymptotics in the limit t ~ * 0, when the quantum sta- 
tistical effects are irrelevant due to low mode occupancies. 
The asymptotic values (dashed lines) are found from Eqs. 
(jTUJ) and (fTTjl : Wcon — w s and Wcorr — 1.12w s , where 
w s ~ 1.62/i?TF is the width of the momentum distribu- 
tion po(fc) = \g(k)\ 2 / x 2 = \^/p^j2Ji(kR TF )/k\ 2 of the 
actual source - the molecular BEC. 

As an alternative measure of the strength of atom- 
atom correlations Q, we also calculate the variance of 
number-difference fluctuations for atoms in different spin 
states and with equal but opposite momenta ±fco, 



V ko ,-k (t) = ([A(N 1M - N 2 .^ ko )} 2 )/SN, 



(12) 



where SN is the shot-noise level that originates from 
uncorrelated states. The number operators are de- 
fined by Nj t ±k (t) — J K dkrij(k,t) [with rij(k,t) = 



Figure 2: (a) Width of the BB and CL correlations relative to 
the momentum width of the molecular BEC, w s ~ 1.62/i?TF, 
as a function of time, for the physical parameters of Fig. [TJ 
(b) Relative number variance Vk a ,-k a (t) as a function of time, 
for R$$, = 250 ftax (1), R^l = 167 ^m (2), and = 83 fim 
(3). The counting length is Afc = tt/2R>^ in all cases. 



(ij(k,t)aj(k,t)], where K is the counting length around 
±fc . On a computational lattice the simplest choice 
that does not require explicit binning of the signal is 
K = Afc, where Afc is the lattice spacing, and therefore 
Nj t ±k (t) — rij(±ko,t)Ak. We emphasize that SN is dif- 
ferent for bosons and fermions. For the bosonic case, SN 
is given by the sum of variances of the individual mode 
occupancies with Poissonian statistics (as in the coher- 
ent state), implying that SN = {Ni l j eo )+ (^2,-fc )- F° r 
the fermionic case, the sum of individual variances gives 
SN = (N hk0 )(l - (N 1)ko ))+ (iV 2 ,_ fe0 )(l - (N 2 ,-ko)) S- 
The variance (fL2l) can be rewritten as 



Vk a ,-k (*) 



1 



Afcni(fcn, t) 



1 — sAkni(ko, t) 
x[g^(k ,-k ,t) - gW(k ,k ,t)},(13) 



where s = 0(1) for bosons (fermions), and we have taken 

^ ^ (2) 

into account that {Ni :ko } = (iVa, fco / an d 9ii v°o> ^o, = 

(2) 

9 22 (~ feoi — &0) t). Variance V ko ,- ko (t) < 1 implies squeez- 
ing of fluctuations below the shot-noise level and corre- 
sponds to a violation of the classical Cauchy-Schwartz 
inequality with g^(ko, —k ,t) > g^ (ko, ko,t) 

The short time asymptotics for the variance can be 
found using Eqs. (fL0|) and (|lip . yielding 



V ko ,- ko (t « t ) = 1 - 37rAfci? TF /32. 



(14) 



The small geometric prefactor in the second term, to- 
gether with the resolution requirement Ak < 1/Rtf, en- 
sures that V ko - ka > 0. We see that the squeezing is 
stronger for larger condensates and counting lengths. 

In Fig. EIL) we plot the variance V kch - ko (t) for three 
different sizes of the molecular condensate. The solid 
lines are the numerical results from Eqs. (jJJ, whereas the 
horizontal dashed lines are the short time asymptotic re- 
sults of Eq. lfl4|) matching precisely the numerical results 
in the limit t — > 0. As we can see, the squeezing of the 
relative number fluctuations (from the numerical curves) 
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does not change significantly with time for bosons, while 
for fermions its dynamics is affected by a stronger depen- 
dence of the fermionic shot noise on the mode occupancy. 
The squeezing is stronger for larger condensates, but is 
still far from perfect squeezing, Vk ,-k (t) = 0, which 
follows from the idealized uniform models @, 0] • 

Gaussian density profile. - For comparison, we also 
give the analytic results for a Gaussian density profile 
po(x) — po exp(— x 2 /2a x ) of the molecular BEC, giving 
the momentum distribution of po(k) oc exp(— k 2 /2af.), 
where a x and <Tk = l/2a x are the rms widths. In this 
case, the atom-atom correlations are 

g$(k,k',t « t ) ~ 1 + 2 e -( fe + fe '> 2 /^/ (< Vpo), (15) 



M. K. Olsen, and C. M. Savage for useful discussions. 



(2) 



(fc,/fc',t«i )-l±e- (fc - fc ' )2/4 ^. (16) 



The respective rms widths are given by dcorr = &k and 
0con^ = V2ak, resulting in the ratio Ccorr /o^rr = V2- 
The relative number squeezing is Vk -k (t <C to) = 1 — 
^/2/-7r Afc a x , with Afc < l/2a x ensuring Vk„,-k > 0. 

Results in 3D. - For a nonisotropic TF parabolic den- 
sity profile of the molecular BEC, performing the inte- 
grals in 3D versions of Eqs. (J6J and (7J) is more cum- 
bersome and we only give the final results for corre- 
lations corresponding to the displacement between the 
pairs of momenta along one of the Cartesian coordinates 
a = x,y, z. The BB and CL correlation widths obtained 

and Wcoir'a — 1.08w s , a , where 
The relative number variance is 
- 15(5* i?^) 3 /2 10 , 



l (BB) 

here are Wcorr a 



w s , a 

Vko,- 



Ws.a 
~ 1.99/i?TF,c 
ko (t < t ) Si 1 



where Afc = 

(Ak x Ak y Ak z )^ 3 and 7^ = (^tf^tf^tf^) 173 are 
the geometric means. The much smaller (than in ID) 
geometric prefactor in the second term, together with 
Ak a < l/i?TF,a, explains why the "raw" (unbinned) 
squeezing is much weaker in 3D @, S| than in ID. There- 
fore, the prescription of Ref. [§] to perform binning for 
obtaining stronger squeezing is more crucial in 3D. 

For a Gaussian density profile, the generalization to 
3D is straightforward. In particular, the BB and CL cor- 
relation widths for a displacement along x, y, or z are as 
in ID, while the relative number variance is Vk ,-k (t ^ 
to) = 1 - \\[2pK Afco 7 ] 3 , where a = {<J x a y a z ) 1 l' i . 

In summary, we have studied the dissociation of a BEC 
of molecular dimers into correlated fermionic and bosonic 
atom pairs. We have obtained explicit analytic results for 
the width and strength of the atom-atom correlations and 
for the relative number squeezing in the short time limit, 
using realistic density profiles of the molecular BEC. The 
results show how the squeezing improves with the larger 
size of the molecular condensate, and how it can degrade 
in strongly inhomogeneous systems. Our approach can 
be easily generalized to describe similar effects in atomic 
four- wave mixing via BEC collisions [2, H, H3| • 
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